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Abstract 

At low temperature the low end of the QCD Dirac spectrum is well described by chiral random 
matrix theory. In contrast, at high temperature there is no similar statistical description of the 
spectrum. We show that at high temperature the lowest part of the spectrum consists of a band 
of statistically uncorrelated eigenvalues obeying essentially Poisson statistics and the correspond- 
ing eigenvectors are extremely localized. Going up in the spectrum the spectral density rapidly 
increases and the eigenvectors become more and more delocalized. At the same time the spectral 
statistics gradually crosses over to the bulk statistics expected from the corresponding random ma- 
trix ensemble. This phenomenon is reminiscent of Anderson localization in disordered conductors. 
Our findings are based on staggered Dirac spectra in quenched SU(2) lattice simulations. 
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The spectrum of the QCD Dirac operator contains important information regarding the 
properties of strongly interacting physical systems. Statistical properties of the spectrum 
completely determine the bulk thermodynamical observables. The low end of the Dirac 
spectrum is particularly important since that part dominates the quark propagators. The 
statistics of low Dirac eigenvalues is fundamentally different in the low temperature chirally 
broken phase and the high temperature chirally symmetric phase. According to the Banks- 
Casher formula [lj the spectral density around zero is proportional to the order parameter 
of chiral symmetry breaking. Therefore in the chirally symmetric phase it vanishes whereas 
in the broken phase it is non-zero. In the latter case random matrix theory (RMT) provides 
an essentially complete statistical description of the low lying part of the Dirac spectrum 
In the intermediate volume, so called epsilon regime the RMT behavior is well understood 
analytically through an effective sigma model description and numerically through lattice 
simulations. 

In contrast, above the finite temperature transition, T c , in the chirally symmetric phase 
there is no well understood statistical description of the low eigenvalues of the Dirac oper- 
ator. In this regime, lacking any analytical insight, one can regard the Dirac operator as a 
randomly fluctuating matrix of size going to infinity in the thermodynamic limit. From this 
perspective there are two fundamentally different possibilities for the spectrum of the Dirac 
operator. If typical random fluctuations can freely mix eigenvectors, eigenvectors become 
extended and the eigenvalue statistics is described by the corresponding RMT. If on the 
other hand, fluctuations in the matrix elements cannot mix the eigenvectors, in some basis 
they become localized and the eigenvalues become independent, obeying essentially Poisson 
statistics. Lattice simulations can test which scenario happens in reality. 

Above T c the spectral density vanishes at zero and RMT has predictions for the eigenvalue 
statistics at such a "soft edge" [3J. Lattice simulations, however, did not find agreement 
with these predictions 0, 5]. Another earlier study 6] did not focus on the spectrum edge, 
but considered full Dirac spectra and found bulk RMT statistics for the full spectrum. 
More recently the possibility of Poisson eigenvalue statistics was suggested again in Ref. 
7|. Based on lattice simulations the authors argued that around T c the low temperature 
RMT statistics is gradually deformed towards Poisson statistics. Ref. [8|, based again on 
lattice simulations, found that although low lying Dirac eigenvectors become very localized 
above T c , the localization is not stable and most likely will disappear in the thermodynamic 
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limit. This finding would disfavor the appearance of Poisson statistics in the spectrum. Very 
recently, using overlap fermion lattice simulations, Ref. |o| found Poisson behavior for the 
lowest two eigenvalues. 

In the present paper we offer a complete understanding of this rather unclear situation. 
We show that above T c the lowest part of the spectrum consists of statistically independent 
eigenvalues obeying Poisson statistics and the corresponding eigenvectors are extremely lo- 
calized. Going up in the spectrum the spectral density rapidly increases and the eigenvectors 
get delocalized. At the same time the spectral statistics gradually crosses over to the bulk 
statistics expected in the corresponding random matrix ensemble. We also show that the 
number of Poisson type eigenvalues depends only on the physical temperature and the phys- 
ical spatial volume and not on the lattice spacing. 

The phenomenon we report here is analogous to Anderson localization occurring in crys- 
talline conductors in the presence of disorder. In that case disorder causes the appearance of 
localized electron states at the band edg e, but for sufficiently weak disorder states towards 
the band center remain delocalized 10]. The corresponding eigenvalue statistics changes 
rom Poisson around the band edge to random matrix statistics towards the band center 



111 ] . It is interesting to note that while in disordered conductors localization is most often 



understood in terms of disorder in the diagonal (on-site) matrix elements of the Hamiltonian, 
in the case of QCD, disorder is entirely in the hopping terms. In fact the on-site matrix 
elements of the staggered Dirac operator are identically zero and all "disorder" is in the 
gauge field in the hopping terms. 

The low lying QCD Dirac spectrum is known to depend strongly on the temporal 
fermionic boundary condition which is effectively a combination of the Polyakov loop and the 



explicitly chosen anti-periodic boundary condition [12] . In the quenched SU(2) theory the 
Polyakov loop Z(2) symmetry is spontaneously broken above T c . Although in the quenched 
theory the two sectors are equivalent, here we only use configurations in the "physical" 
Polyakov loop sector, the one that would survive in the presence of dynamical fermions. We 
use fermion boundary conditions that are anti-periodic in the time direction and periodic in 
all the spatial directions. 

At first we summarize the details of the numerical simulations. The data is based on 
quenched simulations of the SU{2) gauge theory with Wilson plaquette coupling j3 = 2.6 and 
time extension N t = 4. This corresponds to a temperature of T = 2.6T C , well above the finite 
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temperature phase transition. The simulations were done at four different spatial volumes, 
Ng = 16 3 , 24 3 , 32 3 , 48 3 . To assess what happens in the continuum limit we included an 
additional N t = 6, N s = 36, /3 = 2.725 ensemble matched to the N s = 24 simulation in terms 
of physical box size and temperature, but on a 1.5 times finer lattice. On these ensembles we 
computed the 256 smallest positive eigenvalues of the staggered Dirac operator. Due to the 
exact twofold degeneracy of the eigenvalues this yields 128 independent positive eigenvalues 
per configuration that we used for the the statistical analysis. 

As a first step we studied the spatial localization of these eigenmodes. A possible way to 
measure that is through the participation ratio 



where ifi is the normalized eigenvector. If the eigenvector spreads uniformly in a four- 
volume V and is zero elsewhere then its participation ratio is V. Assuming that at this high 
temperature the low eigenvectors can maximally spread in the time direction one can define 
a length scale 



characterizing the spatial extension of the eigenvectors. In Fig. [T] we plot how this quantity 
changes in the spectrum. The averages were calculated in non-overlapping spectral windows 
separately for different spatial box sizes. It is apparent that the lowest eigenmodes up to 
about Xa < 0.22 are very localized and their spatial extension is independent of the box size. 
Above that point the eigenvectors rapidly start to delocalize and their spatial size becomes 
dependent on the box size. This is the point in the spectrum that is known as the mobility 
edge in the context of Anderson localization. 

Besides the localization of eigenvectors the other important factor determining how easily 
eigenvectors can mix is the spectral density. A useful quantity that reflects the combined 
effect of localization and spectral density can be defined as follows. For each eigenvector 
the participation ratio is an approximate measure of the four-volume occupied by the given 
eigenvector. We call the cumulative volume fill fraction the sum of the volumes occupied 
by all the eigenvectors corresponding to eigenvalues less than A. It is understood to be 
normalized by the total box four-volume. In Fig. [2] we show the volume fill fraction as a 
function of the eigenvalue. It is a rapidly increasing function of A. In the lowest part of the 
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FIG. 1: The average linear extension of eigenvectors as a function of the corresponding eigenvalues 
in different spectral windows. The linear size is based on the participation ratio. The different 
curves correspond to different spatial box sizes N s = 16, 24, 32, 48. 

spectrum where the spectral density is small and eigenvectors are localized the fill fraction 
is much smaller than unity. These eigenvectors have very little spatial overlap and they are 
essentially produced independently in different sub- volumes. The corresponding eigenvalues 
are expected to be independently distributed and obey Poisson statistics. In contrast, above 
Xa > 0.28 the volume fill fraction is much bigger than unity and the eigenvectors here 
strongly overlap. Therefore they can freely mix and from this point up in the spectrum the 
eigenvalue statistics is expected to be described by random matrix theory. Between the two 
extremes there must be some transition from Poisson to random matrix statistics. 

A simple way of testing these expectations is to consider the unfolded level spacing 
distribution computed in spectral windows located in the three above described regimes. 
Unfolding is a simple mapping of the eigenvalues that eliminates all the information about 
the spectral density which is not universal, but keeps universal eigenvalue correlations intact. 
Numerically we unfolded by ordering all the eigenvalues by magnitude on all configurations in 
the given ensemble and mapping each eigenvalue to its rank order normalized by the number 
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FIG. 2: The cumulative volume fill fraction for the 24 3 x 4 ensemble. The unfolded level spacing 
distribution will be computed separately in the four shaded spectral windows marked by (a)-(d) 
(see Fig. E]). 

of configurations in the given ensemble. By construction the eigenvalues transformed in this 
way have constant unit spectral density. 

If the original eigenvalues are uncorrelated, the unfolded level spacing distribution is ex- 
pected to be a simple exponential P(x) = e~ x . If, on the other hand, the original eigenvalues 
obey random matrix statistics the unfolded level spacing distribution should follow the so 
called Wigner surmise of the corresponding random matrix ensemble |2(. The staggered 
Dirac operator with fermions in the SU(2) fundamental representation correspond to the 
chiral symplectic random matrix ensemble and the Wigner surmise in that case is [l3| 

= ex P (~f^ 2 ) • (3) 

In Fig. |3] we plot the unfolded level spacing distribution averaged separately for the four 
spectral windows indicated in Fig. |2j For comparison we also show the exponential and 
Wigner surmise distributions expected if the level statistics is Poisson and random matrix 
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FIG. 3: The panels show the unfolded level spacing distribution in different regions of the spectrum. 
The labeling (a)-(d) corresponds to the spectral windows indicated in Fig. [2] with the shaded areas. 
The curved lines are the exponential distribution and the Wigner surmise. 
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respectively. Going upwards in the spectrum the transition from Poisson to random matrix 
statistics is obvious. 

We note that the authors of Ref. also computed the unfolded level spacing statistics 
between the first and the second, the second and the third, etc. eigenvalues separately and 
found that for the lowest few eigenvalues there was a slight tilting of the distribution to the 
left compared to the Wigner surmise (like the one seen in our Fig. |3b). However, they could 
not see the Poisson statistics we report here because they used much smaller spatial volumes 
than what we are using here. Since the spectral density in the Poisson regime is small, in 
a small volume even the first few eigenvalues are in the random matrix regime most of the 
time. 

So far we discussed level spacing statistics for a given spatial box size (N s = 24) and a 
given coarseness of the lattice. For these results to represent real physics it is important to 
check what happens in the thermodynamic and in the continuum limit. In the remainder of 
the paper we address these two questions. 

Thermodynamic limit — We have already seen that the lowest eigenmodes are very lo- 
calized and their spatial extension is not affected by the volume of the box. Therefore we 
expect that these modes occur independently and their average number is proportional to 
the spatial volume. More generally eigenvalues following any intermediate statistics between 
Poisson and random matrix are also expected to occur in numbers proportional to the spa- 
tial volume. If this is true the statistics of eigenmodes number ri\ through m\ in a spatial 
volume V\ should be the same as that of eigenvalues ni through ni2 in a spatial volume 
V2 provided that ^ = ^ = ||. We verified this by comparing the unfolded level spacing 
distribution in two different spatial box sizes, N s = 24 and 32. In Fig. we compare the un- 
folded level spacing distribution of eigenvalues 10-20 in the smaller volume and 24-47 in the 
bigger volume and find good agreement. To demonstrate that these eigenvalues fall in the 
intermediate statistics regime we also indicated in the Figure the exponential and Wigner 
surmise distributions with thin solid lines. This shows that the density of eigenvalues with 
a given intermediate statistics scales with the spatial volume as expected. 

Continuum limit — To see how this picture changes towards the continuum limit we con- 
sidered an additional ensemble on a finer lattice. The parameters of that, N t = 6, N s = 
36, /5 = 2.725 were tuned to match the physical temperature and spatial volume of the 
24 3 x 4 ensemble. Since the average smallest eigenvalue changes roughly with the smallest 
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FIG. 4: The unfolded level spacing distribution for two different spatial box sizes at fixed Nt (a) 
and for two different NfS at fixed physical temperature and physical box size (b). 

Matsubara frequency, it does not make sense to match the spectral windows of these ensem- 
bles. Instead we computed the unfolded level spacing distribution of eigenvalues 10-20 on 
each configuration for both ensembles. As can be seen in Fig. Hb there is perfect agreement 
between the two ensembles. 

On the one hand, this means that the deformation of the distribution from Poisson 
to random matrix statistics occurs through a universal path independently of the lattice 
spacing. On the other hand, it also implies that the number of very localized Poisson 
eigenmodes is proportional to the physical spatial volume and not the volume in lattice units. 
It could suggest that these eigenmodes are localized on some physical gauge field objects that 
survive the continuum limit with a finite physical density. This scenario is also supported 
by the fact that the localization range of the eigenmodes in physical units, computed from 
the participation ratio was roughly the same on the coarser and finer ensemble. 

Discussion — Mixing instanton anti-instanton zero modes could be a natural candidate to 
explain the localized modes. Using the number of zero modes of the overlap Dirac operator 
on these ensembles we estimated the density of topological objects. We found that their 
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density is more than an order of magnitude smaller than what would be needed to explain 
the Poisson modes. No matter what their origin is, however, it is clear that the fermion 
boundary condition plays a crucial role in their appearance because in the opposite SU(2) 
Polyakov loop sector these modes are completely absent |l4|. We hope to return to a more 
detailed discussion of the role of the boundary condition in a later publication. 

Another important issue is how universal the appearance of Poisson modes is. In Ref. {^J 
in a range of spatial volumes the lowest two eigenvalues of the overlap Dirac operator were 
seen to be Poissonian. In a much smaller statistics study than reported here with staggered 
fermions we also verified that the overlap spectrum has similar behavior of crossing from 
Poisson to orthogonal random matrix statistics, the RMT corresponding to the symmetries 
of the overlap operator. It means that the phenomenon we found does not rely on a particular 
lattice fermion formulation. 
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